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AN EXISTENCE THEOREM FOR PERIODIC SOLUTIONS* 

BY 

WILLIAM DUNCAN MACMILLAN 

Introduction. 

Under certain circumstances differential equations which are analytic 
functions of a parameter and of the dependent variables admit of periodic 
solutions which are expansible as converging power series in the parameter. 
The proof of the existence of the periodic solutions is usually much more 
laborious than the actual construction of their first terms. Heretofore it has 
been necessary to give the proof of their existence in each specific case before 
the validity of the construction was assured. 

In this paper the problem of proving the existence of the periodic solutions 
in a very general class of cases will be developed in such a manner that it is 
closely related to that of their practical construction; and it will be shown 
from these relations that, whatever the multiplicity of the solutions, when 
the formal construction is possible the existence proof can be made. Therefore 
in this class of cases the direct existence proof is no longer necessary, and the 
complete discussion is correspondingly simplified and shortened. 

The Differential Equations. 
Let us consider the differential equations 

(i) -J ? = /»(2/i J ■•■, y n , m; (* = i, •••,«), 

where the /,• are uniform, continuous, and periodic functions of t . Without 
loss of generality the period can be taken equal to 2x. Now suppose that 
when the parameter ix is zero, equations (1) admit the periodic solution 

yi=y ( P(t) «-i,".,n), 

in which the period is 2x, or a multiple of 2x. If the period is a multiple of 
2x, we may change t by such a constant factor that it reduces to 2x . 
In order to discuss the solutions of (1) when ju is distinct from zero, let 

* Presented to the Society at Chicago, April 29, 1911. 
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yi = yV + Xi (i-i,-..,n). 

Suppose the /* are analytic functions of ?/,- — t/j 0) = a;, and m > and regular at 

vt - y? = m = o 

for < < ^ 2t . We notice that the right members of the differential equa- 
tions resulting from this transformation vanish for 

xi = ft, = . 

If for convenience of notation we denote the parameter n by Xo , the right 
members of equations (1) can be expanded as power series in the Xj with 
coefficients which are periodic in t, as follows: 

J„ n n j n j k 

(2) f-*+2><*\= ft+ZI^+ISI^w+- 

»* J=l j=0*=0 J=0 teO 1=0 

a = i, •••,«), 

where the 0^?.. are continuous periodic functions of t. 

Before integrating equations (2) let us consider the equations 

(3) d H+£/;\. = o, 

which are obtained by setting the left members of (2) equal to zero. These equa- 
tions are linear and homogeneous with periodic coefficients. A single equation 
of the nth order with periodic coefficients was treated by Floquet,* and the 
set (3) was discussed in detail by Motjlton and MacMillan.I The solu- 
tions of these equations can be written in the form 

(4) Xi=T,A$ ) (t) «-l,-,n), 

1=1 

in which the Aj are constants of integration, while the £p can be so taken that 
the determinant 

A= Uf (0)| 

equals unity. The determinant A expanded according to the elements of the 
ith line is 

(5) A=ZA«|f(0), 

j=i 

where A<*' is the minor obtained by suppressing the ith line and jth column of A 
and attaching the proper sign. 

*Annales de l'Ecole Normale SupSrieure, vol. 12 (1883-1884), p. 47. 
t American Journal of Mathematics, vol. 33 (1911), pp. 63-96. 
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Integration of the Differential Equations as Power Series in the Initial Values. 

If the initial values of the a?» be denoted by at , and if xo be denoted by ao , 
it follows from Cauchy's existence theorem* that the solutions of (2) are 
expansible as power series in the aj (j = 0, • • -, n) , reducing for t = to 

Xi(0) = en; 

and it follows from PoincarS's extension f of Cauchy's theorem that these 
solutions converge for all values of t in any preassigned range for which the 
expansions (2) converge, provided the moduli of the aj are sufficiently small. 
That is, the solutions of (2) can be written 

n 

(6) Xi = 53 x T a j + higher degree terms ( i = 1 , ••• , n), 

0=0 

where 

xf(0) = 1 (i-1, •••, n), 

and all the remaining coefficients vanish at t = . 

On substituting (6) in (2) and rearranging as power series in a , • • • , <x n , 
we have from the coefficients of a } - 

(7) ^ ) +E*i i) 4* )=s0 (i. J-l. ••-. n). 

at k=i 

These equations are the same as (3), and the solutions are therefore 

(8) xf = £,A<]pti ) (i) (i,i-l, -,n). 

From the initial conditions it follows that 

xf(0)= 8 {J , 
where 

5,7=0 (i+i), 5,7=1 (*-j). 
Hence 

Since the determinant of the left members, A , is unity, the solutions of these 

equations are 

Ai = A</\ 
so that 

(9) xf = XX'W) (i, i = 1, ■•• - n). 



* Collected, Works, 1st series, vol. 7. 

ties Methodes.Nouvelles de la Micanique Celeste, vol. 1, p. 55. 
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Since xo = ao, the coefficients of a<> are linear but not homogeneous. From 
equations (2) we find 

ao) -#+£« = ceo c*-i, -.»). 



eft 



*=i 



The general solution of these non-homogeneous equations can be found by 
the method of variation of parameters. We can denote it by 

(11) a£> = m, (t) + t,A«W (* = 1, ••• , n), 

where the m,-(f) are particular solutions, and the ^4 C 4 0) are constants of 
integration. 

From the initial conditions we must have 

4°(0) = o. 

Hence 

(12) iU?W) = -!M0) (i = i, ••-,«). 

Solving these equations, we get 

4 0) =-tAfm,-(0) (fc-1. ••-.n), 

so that 
(13) xf=m i (t)-t > t,^m j (0)^(t) (*-l. ■■■.n). 

The sum of all of the linear terms in the solution for the Xi is found from (9) 
and (13) to be 

£ *?<*! = E E W kf aj-iti, AT m, ( ) tf> a + m, ( < ) • « 
(14) 

= ZrEA« ) [a,-m,.(0)a ]l|i i >(<) + m j (0-ao- 

A Change of Parameters. 
Let us now make a change of parameters by the linear substitution 

(15) «o = jSo, J2^[a i -m i (0)a ] = Ph (* = 1. ••-, n), 

which can be solved for the as in terms of the j3's, since the determinant 
| A^' | has the value unity. The result is 

n 

(16) «o = 0o, ctj = my(0)|8o + £#>(<)) 0* (i - 1. ••-, n). 

Trans. Am. Matb. Soc. 10 
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Since the solutions (6) as power series in a , •••,«» converge provided the 
moduli of the a,- are sufficiently small, they can also be developed as power 
series in the ft which will converge provided the moduli of the /3»- are suffi- 
ciently small; and for t = will reduce to 

n 

x,(0) = m,-(0)j8o + E6°(0) ft «-l. ••■.«)• 

Since the transformation (15) is linear and homogeneous, the coefficient of 
every /? of degree higher than the first will vanish at t = . Expressed in 
terms of the fts, the solutions can be written 

*< = mt (0 /?o + £ #> (0 ft + £ E^ftft 

(17) n ^ 4 

,,=o*=o «=o 

Since ft , • • • , ft, can be regarded as arbitrary, one sees that they are nothing 
else than the constants of integration of the linear terms (3). We have 
therefore, since the periodic properties of (1) have not been used, the 
Theorem: If 

dx ■ 

(A) -^=fi(x 1 ,---,x„;n;t) (* = 1, •••,») 

is a system of differential equations in which the ft are analytic functions of 
%\, • • • , %n', M> regular at the point Xi = • • • = x n = fi = 0, for all t in the 
interval ^ t < T , and vanish for X\ = • • • = x n — n = , and if the fi are 
uniform and continuous in tin the interval ^ t ^ T; and if 

(B) f + I>^-=0 «-i. ...,») 

are the linear terms of (A) equated to zero, then the solutions of (A) are expansible 
as power series in fi and the constants of integration of the solutions of (B), and 
these solutions converge in the interval ^ t ^ T provided the moduli of the 
parameter n and the constants of integration are sufficiently small. 
For the sake of future notation we will write (17) in the form 

(is) x t = £*?%+££ x%^ k + £ £ z x% ^.^ l + ■■-, 

though it is to be observed that the x^p are the same as the linear terms of (17), 
viz: 

4 i)= #' («.i = i. •■•■ «). xf=m i (* = 1, ••-,«), 

and are not the same as the a^* 5 used in (6)-(14). 
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On substituting (18) in (2) and rearranging as power series in (3 , • • • ,Pn, 
we find from the coefficients of the terms of the second degree 

j=0*=0L aJ *=1 J j'=0*=0 \r=0 / \«=0 / 

(19) = ± ± 6„ (£ £ W^+xfxf)) PA 

r=0 3=9 \j=0 k=0 / 

= He J H C(4X' } +4'H r) )) **A> 

,7=0 *=0 \r=0«=0 / 

where e,* = 1 if j 4= A; , and e# = 1/2 . The terms 

r=0 8=0 

are all known as soon as the linear terms are known. Their sum depends only 
on the indices i, j , k. Consequently the particular solution of (19), which 
depends on this sum, can be denoted by ip%{t) , and the solution of (19) is 

(20) ss^ft'iiU^o+s^'kft. 

j=0*=0 j=0J=0L !=1 J 

By virtue of the initial conditions we must have 

Z4^°(0) — *$(o) «-i. •••-«). 

whence 

so that (20) becomes 

(2D £ £ JMP k = £ £ [^(o - £ £ a?'^(o)^(o1 <w 

j=0 4=0 J=0 *=0 L i=l A»=l J 

In the same manner we obtain for the terms of the third degree 

III -f+I« W, 

w p q / n \ / n \ / n \ n p 

= ziie z^- i4^ i^ +^^^ 

j>=0 J=0 »-=0 \ 3=0 / \ k=0 / \l=0 / p=0 g=0 

l \ j=0 / \ j=0 *=0 / \ .7=0 / \ j=9 *=o / J 
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n j k 

j=0 k-.Q 1=0 



n p q 

2Ll 2^ Zw *jkl Vpqr 
p=Q q=q r=0 



+ ZI»,««B' 



p=0 q=0 



PifoPl 



+ (x ( J " ) x ( * ) +x\ p) x<f ) )x'-; ) 
+ {x^xf+^xf)*p 

+ (x$xp+xf k xr) 

+ (x^xf+x ( f>x^) 
+ {x<$xf + xf l x?>) 

The symbol tju means 1 ii j , k and I are all different; 1/2 if two are alike and 
different from the third; and 1/6 if all three are alike. The symbol 5^2 means 
that when two of the three letters j , k , and I are equal only one of two like 
terms is to be taken, and when all three are equal only one of three like terms 
is to be taken. 

The right members of equations (22) are all known and they depend 
only upon the indices i, j , k, and I . If therefore we denote the particular 
solutions by <?$,(<), the complete solutions can be written 

(23) e e i: 4'iftftis, = e e £ [ *$ a ) + i, a% &> I w k & . 

j=0 4=0 1=0 j=0 *=0 1=0 L m=l J 

On determining the constants of integration A^ k \ so as to make the a^\ vanish 
at t = , we get 

n j k n j k l~ n n ~1 

(24) SII x% Wl = E E E ]<?% ( o - E E Af ^ ( o ) #> ft./s^, 

j=0WH ,-=Jfe0i=nL ST=1A=1 J 

(i = 1, •••, n). 

The coefficients of the fourth and higher degree terms in /3i , • • • , /3„ can be 
developed in a similar manner. 

Integration of the Differential Equations as Power Series in /*. 

Without specifying the initial values of the xi and without regard to the 
convergence of the series so derived, the differential equations can be integrated 
formally as power series in fi , or any root of ju , say fi llp . We will integrate as 
power series in fi, and leave the constants of integration arising at each step 
undetermined. We shall have then 



(25) 



x t = atfV + 4V + 4V + 



(i = l. 



n). 



Substituting these expressions in (2), we have from the coefficients of the first 
power of ii, since x = n , 



(26) 



dxf 



dt 



h+E<W=^(0 



(*-l, 



, n). 



3=1 



(i = 1, ■■■,n). 
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These equations are the same as equations (10). They have therefore the 
same solution (11), viz., 

(27) x[ l) = m, (0 + Z 4 ; W } ( * = l ' • • • • n >- 
which can be written 

(28) xf = £ ^a-f (» = 1 , ■ • • . n) , 

,?=o 

the x ( p being the same functions of t as in (18), and A^ = 1 . 
The coefficients of ju 2 give the equations 

(29) ^ 2 - + z w =zi^ z Af^ ) ( z « ) 

These equations are essentially the same as (19), the right members differing 
only in that ^° is substituted for /?, . The solutions are therefore 

(30) xf=±i:< P f k {t)AfA?+±B?tP « = i. -.«)• 

j=0 4=0 J=l 

We propose to leave the constants of integration undetermined, but we observe 
that if we give the Bf } the form 

Bf= -.±A« ±£<pfl(0)A?A?+Af>, 

h=l j=0 *=0 

where the ^4 ( , 2) are undetermined, the solution (30) takes the form 

xf = z z \ v % a) - £ £ a?v?2 (o) #> (o~Lw+ i^fr 

;=0 *=0 L Z=l A=l J J'=l 

(3D 

— £^, £^, Xj k Jl, A h -f- £^, -Clj X, , 

j=0 k=0 j=0 

where A ( 2) = . The initial value of a4° is then 

j=0 

which is undetermined since all of the A ( f are undetermined. 
From the coefficients of /* 3 we have 

J (i) » n p q / n \ / n \ / n \ 

% + Z W = Z Z Z «& ( Z 4'M P) ) ( X ^«> ) ( Z Afxf ) 

«t ,/=l j=0}=0 r=0 \j=0 / \*=0 / \l=0 / 

(32) + Z Z <% { £ Afxf \ £ Z AVAfx% + £ Afxf 1 

p =o «=o I j=o L *=o «=o t=o J 

+ Z ^f [ Z Z AfA<M}+ £ Afx^A \ , 

j=0 L *=0 1=0 *=0 J J 
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which becomes on rearranging the right member 



[April 



J (t) n n j k 



2w 2-1 2—1 e jia"q t 



p=0 q=0 r=0 



jkl qpr 



+ Wfxf+xfW)x? 

+ (x ( / ) xf+xfx\ p) )x^ 

l+ix^f+xfx^x?} 



(33) 



n p 
p=0 y=0 



+ IIy^ 



+ (4W ) +^M P) ) 
+ (4?4 ?) +*^ ( / ) ) 



4 1) 4 1> 4 1) 



+ i: i:[i: t^o4^ ) +^ ( / ) )l[^ , ^ 2) +4 2) ^ l) ]- 

^=o *=o Lp=o ?=o J 

On comparing these equations with (22) and (19), it is seen that the solution is 

n j k n j n 

(34) *?> =ZIE <P%AfA^Af + 11 «$[44 2) + 4 W 4 1) 1 + 22 W • 

,7=0 4=0 J=0 j=0 *=0 y=l 

If now the constants of integration are given the form 



Bf = - £ A?> [ £ £ £ <p%( )AfAfAf 

A=l L J'=0 k=0 1=0 

n j "1 

+ IE ^(0)(4"4 2, +if4") + 4*> 

j=0 4=0 J 

equations (34) become 

» j k I~ 71 71 ~| 

,7=0 *=0 1=0 L £T=1 ft=l J 

+ ££\<pfk-££ a ( ;v^(o)^1 [AfAf+ AfAf] + £ #?$>, 

j=0 k=0 \_ g=\ h=l J .7=1 

which are simply 



(35) 



n j k 
,7=0 k=0 1=0 



+ ZZ^ (44 2) + 4 2) 4 I) ) + 22 Afxf, 

}=0 k=0 )=0 



where the Af } are the constants of integration, A^ = , and xfy and x ( p k 
are the same functions of t which occur in (24) and (21). Since these functions 
vanish with t , the initial values of the a4° are 



x?(Q) = 22Afx?(0), 

j=0 
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The coefficients of the higher powers of m can be determined in a similar 
manner. So far as they have been worked out, the solutions as power series 
in n are 

[n ~l I - re i re -| 

£ Afxf U II x%AfAf + Z xfAf m 2 

(36) +\ttt ^AfA^A? +H 42 UfAf + AfA<?) 

L j=0 k=0 1=0 j=0 4=0 

+ T,A?xf] l S+..., 

5=0 J 

where Af =1 and A { p) = 0, (p = 2 , • • • , oo ) and all the other Af are 
undetermined. But equations (36) are precisely what one obtains on sub- 
stituting in (17) 

(37) ft = /i, 0, = AfiL + AfV + Afn 3 + ■■■ U = 1, ••• , n), 
This, of course, should be expected, for if we take 

0<A</>< M, 

where M is a positive constant, the substitution (37) converges for all values 
of m < 1 , and consequently the substitution of (37) in (17) will give the solu- 
tions of the differential equations as power series in /* which converge for values 
of n sufficiently small, and which reduce at t = to 

x i (0) = T,i,AfW;Xo)S. 

But these are exactly the conditions under which (36) were developed. The 
two series are therefore identical. 

Conditions for Periodic Solutions. 

From the hypothesis as to the Bfcl in the coefficients of the differential 
equations, i. e., that they are periodic functions of t with the period 2ir , it 
follows that sufficient conditions that the solutions shall be periodic with the 
period 2kir ( k an integer ) are 

Xt ( 2kir ) = Xi (0) (i = 1, ••• , n). 

If the X, all return to their original values, it is obvious from the differential 
equations that their first derivatives retake their initial values, and therefore 
all higher derivatives do likewise. Consequently under these conditions 

xt(t+ 2kir) = Xi(t) (» = !,•••,»). 
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If we denote the difference between the value of a function at t = 2kir and 
at t = by a dash over the letter representing the function, e. g., 

Xi = x t ( 2kir ) — Xi ( ) , 

the conditions for periodicity as derived from (18) are 

n n j m j k 

(38) = I^ + H xf k jS.p k + IZI x% fyW, + ■ ■ ■ . 

j=0 j-0k=0 j=0k=0 1=0 

As derived from (36), the conditions are 

Y.Afxf U EE^'^ + Eifif L 2 

j=0 J L j=o *=o i=o J 

[ n j k n k 

j=0k=0 1=0 j=0k=0 

+ ±Afxf^+---. 

Since equations (39) can be derived from equations (38) by the substitution 
(37), it is clear that if the undetermined constants A^ p) can be determined so 
as to satisfy (39), then the values of the /?/ as defined in (37) will satisfy the 
conditions (38). That is, the determination of the constants Af* so as to 
make the series (36) periodic is equivalent to a purely formal solution of the 
equations of condition (38). The convergency of the series thus obtained 
is thus reduced to a question in the theory of implicit functions: Under what 
conditions does a purely formal solution of a set of equations of the type (38) 
converge ? 

Let us suppose first that in the coefficient of n p every A^f can be determined. 
This will be possible if the determinant | ^*' | =f= ( i , j = 1 , • • • , n,. ) , and 
the determination is unique. But if this determinant is not zero equations (38) 
can also be solved uniquely, since | af^ \ is its functional determinant. 

The condition that the determinant 

| if | = \if\ (i.i-i, •••.*) 

shall not vanish is equivalent to the condition given by Poincare,* but is some- 
what simpler. The functional determinant, | d\pi/dfa | in his notation, is the 
functional determinant of the a ( obtained from (38) after transforming them 
by equations (15). It can be shown without trouble that the functional 
determinant of the a,- is the product of the two determinants | if ( ) | and 
| |f | , the first of which (see equation (4) et seq.) has the value unity. The 
second can be formed as soon as a fundamental set of solutions of (3) is known. 
* Les MUhodes Nouvelles de la Micanique Celeste, vol. I, p. 83. 
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Let us suppose that the roots of the fundamental equation for the solutions 
of (3) are all distinct. Then 

where the yp { p ( t ) are periodic with the period 2tt . The determinant 

| £« | = | J,f ( ) I II ( e 2kXi " - 1 ) = II ( e 2% " - 1 ) (* an integer) , 

j=l 3=1 



can vanish if, and only if, k\j = mod V — 1 for one or more values of j . 
The same result holds if two or more of the characteristic exponents Ay are 
equal. Consequently periodic solutions with the period 2kir always exist 
if k\j ^ mod V — 1 for all j . In particular, if k = 1 we have the theorem 
that, if none of the characteristic exponents Ay are congruent to zero mod V — 1 , 
then there exists one and only one periodic solution with the period 2t .* 

It is to be observed that the determinant | |y° | may vanish for t = 2kir 
while it does not vanish for t = 2ir . Indeed, if any of the Ay are pure imag- 
inaries and rational, there exist values of k for which the determinant vanishes. 

If the determinant | £$ l) | vanishes, then the constants A { f cannot all be 
determined by the periodicity conditions on the coefficients of fiP , and in gen- 
eral the solutions will no longer be unique. Nevertheless, it may still be 
possible to determine the Af> so as to satisfy the periodicity conditions, and 
the solutions thus obtained will converge provided the determination becomes 
unique at any step; for it is shown in a paper entitled A method for deter- 
mining the solutions of a set of analytic functions in the neighborhood of a 
branch -point t that a formal solution of a system of equations of the type of 
(38) as a power series in integral or fractional powers of n converges provided 
the solution is not multiple identically in n . If the solution is such a multi- 
ple solution, then at no step will the determination of the constants A) p) become 
unique, and the question of the convergence of such coincident branches re- 
mains open. We have then the 

Theorem: // 

dl = ^ ! ^ Xl ' " ' Xn] M ' ^ (i = 1 ' '"' n) 

is a system of differential equations in which the /, are expansible as power series 
in Xi, • • • , x n , and ix , vanishing for Xi = • • • = x n = n = , with coefficients 
which are uniform, continuous and periodic functions of t with the period 2x; 
and if the /,• converge for ^ t^ 2ir when x, < pi, n<r , then the solutions xt(t) are 

* This theorem was given by Poincare, loc. cit., p. 181. 

t MacMillan, Mathematische Annalen. The paper is in type and will appear 
in 1912. 
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expansible as power series in n, or any fractional power of n , which converge for 
all t in the interval 5^ t ^ 2kir provided \ n | is sufficiently small. If the constants 
of integration arising at each step can be determined so as to make the solution 
formally periodic with the period 2kir, then the solution so determined will be 
periodic and converge for all finite values of t provided \n\is sufficiently small. 

If periodic solutions exist as power series in fractional powers of n , but not 
in integral powers, then it will not be possible to determine the constants A { p 
in (36) so as to satisfy the periodicity conditions, but on taking 

Afn = ft (j-l, ■••, n), 

A ( j p) =0 (p = 2, •••<»), 

we pass at once to the general equations of existence (38) in ft , • • • , ft, , which 
are more convenient and much simpler than the corresponding equations 
expressed in terms of the initial values of the Xi . The forms of the periodic 
solutions, if they exist at all as power series in integral or fractional powers of 
n, are then readily found. But if the construction of the solutions can be made 
in either integral or fractional powers of fi, one is assured of the convergence 
of the solutions so obtained for all values of t provided the modulus of n is 
sufficiently small. 

The University op Chicago, 
June 10, 1911. 



